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Abstract 

We show how to transform the problem of finding d + 1 mutually unbiased bases in 
into the one of finding d{d+ 1) vectors in C^^. The transformation formulas admit a solution 
when d is a prime number. 
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1 INTRODUCTION 



The determination of mutually unbiased bases (MUBs) is of pivotal importance in the theory of 
information and in finite quantum mechanics. Let us recall that two orthonormal bases of a unitary 
space are said to be unbiased if the modulus of the inner product of any vector of one basis with 
any vector of the other is independent of the considered vectors (see Section II for a definition of 
MUBs in C"'). Such bases are useful in classical information (network communication protocols) 
im and quantum information (quantum state tomography and quantum cryptography) @ as well 
as for the construction of discrete Wigner functions |[3l, the solution of the mean King problem 
JH and the understanding of the Feynman path integral formalism [5J. They are at the root of a 
formulation of the Bohr complementarity principle for finite quantum systems. 

There exist numerous ways of constructing sets of MUBs. Most of them are based on discrete 
Fourier transform over Galois fields and Galois rings, quadratic discrete Fourier transform of 
qudits, discrete Wigner distribution, generalized Pauli operators, generalized Hadamard matrices, 
mutually orthogonal Latin squares, finite geometry methods, angular momentum theory. Lie-like 
approaches, and phase states associated with a generalized Weyl-Heisenberg algebra (see and 
fTl for a review on the subject). 

The aim of this note is to introduce a transformation that makes possible to replace the search 
of d + 1 MUBs in by the determination of d{d+ 1) vectors in C''^. 



2 MUTUALLY UNBIASED BASES IN 

Two distinct orthonormal bases 

Ba = {\aa) : a = 0,1,- ■ -,(1-1} (1) 

and 

55 = {|6/3):/3 = 0,l,...,d-l} (2) 
(with a 7^ 6) of the d-dimensional Hilbert space C'' (d > 2) are said to be unbiased if 

|(a«|6/3)| = -L, (3) 

where ( | ) denotes the inner product in C^. It is well-known that the maximum number of MUBs 
in is d + 1 and that a complete set ofd + l MUBs exists if d is prime or the power of a prime 
number [[Tl, [[8l, flU. On the other hand, it is not known if it is possible to construct a complete set 
ofd + l MUBs in C'' in the case where d is not the nth power (n G N*) of a prime. However, in 
this case there exists at least 3 MUBs, a well-known result for d = 6. In spite of a great number 
of numerical studies, no more than 3 MUBs were obtained for = 6 [HOl, ([HI, [[HI, [[T3l, [fT4l. 
in agreement with the fact that it is widely believed that only 3 MUBs exist for c? = 6. 
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If we include the a = b case, Eq. (3) leads to 



or equivalently 



I (aa 1 6/3) I = 6a^p6a,b + -^(1 - Sa,b) (4) 

\{aamf = 6a,f,6a,i> + ^{l - 6a,b). (5) 

We note the presence of a modulus in Eqs. dH) and ([5]). This modulus certainly constitutes a 
handicap when performing numerical calculations. 

3 PASSING FROM TO C^' 

The problem of finding a complete set ofd + 1 MUBs in amounts to find d{d + 1) vectors \aa) 
satisfying Eq. ©, where a = 0, 1, . . . , (i and a = 0,1, . . . , d—1 (the indexes of type a refer to the 
bases and, for fixed a, the index a refers to one of the d vectors of the basis corresponding to a). 
By following the approach developed in [fTSl for positive operator valued measures and MUBs, 
we can transform this problem into a (possibly) simpler one (not involving a square modulus like 
in Eq. (|5])). The idea of the transformation is to introduce a projection operator associated with 
the \aa) vector. 

Let us suppose that it is possible to find d + 1 sets Ba (with a = 0,1, ... ,d) of vectors of 
such that Eq. © is satisfied. It is thus possible to construct d{d + 1) projection operators 

Hao = |aa)(aa|, a = 0,l,...,d, a = 0,1, . . . ,d — 1. (6) 

From Eq. ([5]), it is clear that the Haa operators satisfy the trace condition 

Tr (n„„nfe^) = 6o^,p6a,b + ^(1 - ^a,b), (7) 

where the trace is taken over C^. Each operator Uaa can be developed on an orthonormal basis 
{Epq : p, q = 0, 1, . . . , d — 1} of the space of linear operators on (orthonormal with respect to 
the Hilbert-Schmidt inner product). In other words 

d-l 

= ^ ^ Wpg{aa)Epg. (8) 

The Epg operators are generators of the GL{d, C) complex Lie group. Their main properties are 

E^g = Egp, EpqErs = Sg^rEp^s, Ti [El^Ers] = 5p,r5g,s, P, q,r,s E Z/dZ (9) 
and they can be represented by the projectors 

Epg=\p){q\, p,qeZ/dZ. (10) 



2 



The Wpg{aa) expansion coefficients in Eq. ([8]) are complex numbers such that 



Wpq{aa) = Wqp{aa), p,qGZ/dZ, (11) 



where the bar denotes complex conjugation. 
By combining Eqs. ^ and ([8]), we get 

d-l d-1 ^ 

'Y^'^Wpq{aa)Wpq{h(5) = Sa,l3Sa,b + ^(1 -Sa,b)- (12) 
p=0 q=0 

The Uaa operators can be considered as vectors 

w{aa) = {woo{aa) , woi{aa) , . . . , WmmicLO:)) , m = d — 1 (13) 
in the Hilbert space C"'^ of dimension d"^ endowed with the usual inner product 

d-l d-l 

w(aa) ■ w(6/3) = ''^^ Wpq{aa)wpq{b(3) (14) 

p=0 q=0 

(in Eq. (fT3l) . we use the dictionary order for ordering of the components of w(aa)). Equation 
(fT2l) can then be rewritten as 

w{aa) ■ w{b(3) = 6a,p6a,b + 3(1 - ^a,^), (15) 

d 

to be compared with Eq. ([5]). 

The determination of the Uaa operators and, therefore, of the \aa) vectors in C^, is equivalent 
to the determination of the Wpq{aa) components of the w{aa) vectors in C"'^. This yields the 
following. 

Proposition 1. For d > 2, tofindd + 1 MUBs in (if they exist) is equivalent to find d{d+l) 
vectors w(aa;) in C^^ satisfying 

w{aa) ■ w(a/3) = 6a, /3 (16) 

and ^ 

Mv{aa) ■ w{bf3) = - for a 7^ 6, (17) 

where a,b = 0,1, . . . ,d and a = 0,1, . . . ,d — 1. 

Proof. The proof follows from Eqs. (l6l)- (fT5l) . □ 

For a b, Eqs. (fT6l) and (fTT]) show that angle Uaabf^ between any vector w(aa) and any vector 

w{b(3) is 

UJaabl3 = COS'\l/d) (18) 
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and therefore does not depend ona, a,b and /3. 

Proposition 1 can be transcribed in terms of matrices. Let Maa be the Hermitian matrix of 
dimension d whose elements are Wpg{aa) with p,q E Z/dZ. Then, Eq. (fT4l) gives 

w{aa) ■ w{b(3) = Tr (M^^M,,^) . (19) 

Therefore, we have the following proposition. 

Proposition 2. For d > 2, tofindd + 1 MUBs in (if they exist) is equivalent to find d{d + l) 
Hermitian matrices Maa of dimension d satisfying 

Tr (MaaMh^) = 6a,^6a,b + ^(1 - ^a,b), (20) 
where a,b = 0,1, . . . ,d and a = 0,1, . . . ,d — 1. 

Proof. The proof is trivial. □ 
Finally, as a test of the validity of Propositions 1 and 2, we have the following result. 
Proposition 3. For d prime, Eqs. f fT^I ) and dTTI) or Eq. ( |20l ) admit the solution 

Wpq{aa) = l^ei-(p-i)i(d-2-P-ci)a-2a]/d^ a,p,qE Z/dZ (21) 

(Jj 

and 

Wpg{da) = 6p^q6p^a, a,p,qeZ/dZ (22) 

for a = d. 

Proof. The proof is based on the use of Gauss sums lfT6l in connection with ordinary [fTTl and 
quadratic [fTSll discrete Fourier transforms. Indeed, it is sufficient to calculate w{aa) ■ w(6/3) as 
given by (fT4l) with the help of (|2TI) and (|22l) in the cases a = b (for a = 0,1, . . . , d), a ^ b (for 
a,b = 0,1, . . . ,d — 1) and a ^ b (for a = 0,1, . . . , d — 1 and b = d). The main steps are the 
following. 

(i) Case a = b = d: We have 

d-l d-1 

w{da) ■ w{d(3) = ^p,q^p,»^p,l3 = (23) 

p=0 5=0 



(ii) Case a = b = 0,1, . . . ,d — 1: We have 

d-l d-1 

d^ 



w{aa) ■ w(a/3) = ^J^Yl e'^^^^-^^^""'')/'^ = 6^,0. (24) 



p=0 9=0 
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(iii) Case a b (a = 0,1, . . . ,d — 1 and b = d): We have 



p=0 (7=0 p=0 



(iv) Case a 7^ 6 (a, 6 = 0, 1, . . . , ci — 1): We have 

Ci2 



w(aa) ■ w(6/3) = E E e^-(P-^)[('^-2-p-g)(;,-a)+2(a-/3)]/d_ ^26) 



p=0 (j=0 

The double sum in (|26l ) can be factored into the product of two sums. This leads to 

(i-l 



w(aa) ■ w(6/3) = ^ 



gi7r{(a-fe)fc2 + [(d-2)(fe-a)+2(a-/3)]fe}/d 



fc=0 

By introducing the generalized Gauss sums [fT6ll 



(27) 



5(m,w,«;) = ^ e'"("^'+"^)/'", (28) 

fc=0 

(where u, f and w are integers such that u and to are coprime, uw is nonvanishing and uw + f is 
even), we obtain 

^{aa) = j^\S{u,v,w)\\ (29) 

with 

u = a-b, v = -{a-b){d-2) + 2{a-f3), w = d. (30) 
The S{u, V, w) Gauss sum in (|29l)-(l30l) can be calculated from the methods in lfT6l . This yields 

w(aa) ■ w(bp) = (31) 
d 

which completes the proof. □ 

4 CONCLUSION 

As a conclusion, passing from to amounts to replace the square of the modulus of the inner 
product {aa\b/3) in (see Eq. ([5])) by the inner product w{aa) ■ w(6/3) in C^^ (see Eqs. (fT6l ) and 
(flTl)). It is expected that the determination of the d{d + 1) vectors w(aa) satisfying Eqs. (fT6l) 
and (fTTI) (or the d(d + 1) corresponding matrices Maa satisfying Eq. (|20|) ) should be easier than 
the determination of the d{d + 1) vectors \aa) satisfying Eq. ©. In this respect, the absence of a 
modulus in (fTTI) represents an incremental step. 

Now we may ask the question: How to pass from the w(aa) to \aa) vectors? Suppose we 
find d{d + 1) vectors of type (fT3l) satisfying Eqs. (fT6l) and (fTTI) . Then, the Haa operators given by 
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([8]) are known. A matrix realization of each Uaa operator immediately follows from the standard 
matrix realization of the generators of the GL{d, C) group. The eigenvector of the matrix of Uaa 
corresponding to the eigenvalue equal to 1 gives the \aa) vector. 

Of course, the impossibility of finding d{d + 1) vectors w(aa) or d{d + 1) matrices Maa 
would mean that d + 1 MUBs do not exist in when d is not a strictly positive power of a prime. 

Transforming a given problem into another one is always interesting even in the case where 
the new problem does not lead to the solution of the first one. In this vein, the existence problem of 
MUBs in d dimensions was approached from the points of view of finite geometry, Latin squares, 
and Hadamard matrices (see [6 J and references therein) with some interesting developments. We 
hope that the results presented here will stimulate further works, especially a new way to handle 
the d = 6 unsolved problem. 

To close, let us mention that it should be interesting to apply the developments in this paper 
to the concept of weakly MUBs recently introduced for dealing in the Z/dZ x Z/dZ phase space 

m. 
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